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ABSTRACT 

It is well accepted that 'hot Jupiters' and other short-period planets did not form 
in situ, as the temperature in the protoplanetary disc at the radius at which they 
now orbit would have been too high for planet formation to have occurred. These 
planets, instead, form at larger radii and then move into the region in which they 
now orbit. The exact process that leads to the formation of these close-in planets is, 
however, unclear and it seems that there may be more than one mechanism that can 
produce these short-period systems. Dynamical interactions in multiple-planet systems 
can scatter planets into highly eccentric orbits which, if the pericentre is sufficiently 
close to the parent star, can be tidally circularised by tidal interactions between the 
planet and star. Furthermore, systems with distant planetary or stellar companions 
can undergo Kozai cycles which can result in a planet orbiting very close to its parent 
star. However, the most developed model for the origin of short period planets is one in 
which the planet exchanges angular momentum with the surrounding protoplanetary 
disc and spirals in towards the central star. In the case of 'hot Jupiters', the planet 
is expected to open a gap in the disc and migrate in, what is known as, the Type II 
regime. If this is the dominant mechanism for producing 'hot Jupiters' then we would 
expect the currect properties of observed close-in giant planets to be consistent with an 
initial population resulting from Type II migration followed by evolution due to tidal 
interactions with the central star. We consider initial distributions that are consistent 
with Type II migration and find that after tidal evolution, the final distributions can 
be consistent with that observed. Our results suggest that a modest initial pile-up 
at a 0.05 au is required and that the initial eccentricity distribution must peak 
at e 0. We also suggest that if higher-mass close-in exoplanets preferentially have 
higher eccentricities than lower-mass exoplanets, this difference is primordial and is 
not due to subsequent evolution. 

Key words: stars: formation — stars: pre-main-sequence — circumstellar matter — 
planetary systems: protoplanetary discs — planetary systems: formation 



1 INTRODUCTION 

The first extrasolar planet detected around a solar-like star 
was 51 Pegasi b, discovered in 1995 ( Mayor fc Qucl oz 1995|). 
This planet, with a mass about half that of Jupiter and a 
semimajor axis of 0.052 au, was the prototype of a class of 
planets now known as 'hot Jupiters'. These are gas giant 
planets that orbit close (a ^ 0.1 au) to their parent stars. 
There is, however, a general consensus that these planets did 
not form in situ since the temperature in the protoplanetary 
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disc at the radii where these no w orbit would b e too high 
for planet formation to proceed (|BeIl. et al.lll997l ). 

The most developed model for the origin of 'hot 
Jupiters' and other short-period planets is one in which these 
planets exchange angular momentum with the surround- 
ing protoplanetary disc and spiral in towards the centra l 
star (|Goldreich fc Tremainelll980l : iLin fc Papaloizoulll986l '). 
In the case of 'hot Jupiters' the planet is expected to open 
a gap in the disc and migrate in what is known as the Type 
II regime. Lower-mass , close-in pl anets may migrate in the 
gapless Type I regime (|Wardlll997l ). This work will, however, 
only be considering planets that could potentially have un- 
dergone Type II migration. 

Migration is not the only mechanism that can lead to 
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planets orbiting very close to their parent stars. Dynamical 
interactions in multiple planet systems - often referred to 
as planet-planet scattering - can scatter planets into highly 
eccentric orbits which, if the pericentre is sufficiently close 
to the star, can be circularised by tidal interaction s between 
the planet and the host star (iRasio et al.lll99^ . Further- 
more, systems with distant stellar or planetary companions 
on inclined orbits (with respect to the inner planet's or- 
bit) can undergo Kozai cycles which, if then followed by 
tidal evolution, can result in a planet orbiting very close to 
its parent star jKozaill 19621 : lEggleton. Kiseleva fc Hu^l 19981 : 
IWufc Murravll2003l ). 

These different mechanism will result in different dis- 
tributions of close- in planets. Early s imulations of planet- 
planet scattering (jPord. Havlickova fc Rasio 2001) found 
that the likelihood of a massive planet being scattered into 
an orbit that brought it very close to the central star was less 
than 1 %, suggesting that this was unlikely to be the primary 
mechanism f or producing 'hot Ju pitcrs'. Longer term inte- 
gration (Mar zari F. fc Weidenschi lline 20o3) increases this 
to ~ 10 %, while iNaeasawa. Ida fc B cssho (2008) suggest 
that "Kozai migration" l|Wu fc Murray. i2003i ) can further 
incre a se this to ^ 30 %. Recent observations (|Winn et al.l 
I2OO9I : ICoUier Cameron et all l2010l : lOueloz. et all l2010l ) of 
exoplanets with orbits that are retrograde with respect to 
the rotation of the their central star suggests that "Kozai 
migratio n" must play a role in the formation of some close-in 
planets (jTriaud et al.ll2010l ). 

Theoretical modelling (jFabrvckv fc Tremaind |2007| ) 
and analysis of the curre nt sample of close-in exoplanets 
l|Morton fc JohnsonI [ioTH ') suggests that Kozai cycles to- 
gether with tidal friction can produce some, but not most, 
of the close-in planets. In this work we intend to investi- 
gate if initial conditions consistent with Type II migration 
of giant planets can lead - after tidal evolution - to distribu- 
tions that are consistent with those observed. Although it 
is fairly clear that Type II migration is not the only mecha- 
nism capable of producing close-in giant planets, if it is the 
dominant mechanism we would at least expect the current 
distribution to be consistent with one having evolved from 
an initial distribution resulting from Type II migration. The 
paper is structured as follow. Section 2 describes the basic 
model and assumptions, Section 3 dicusses the results and 
in Section 4 we discuss the results and draw conclusions. 



2 BASIC MODEL 

In this paper we use knowledge of Type II migration to 
choose initial semi-major and eccentricity distributions of 
gas giant planets around Sun-like stars and evolve these sys- 
tems by integrationg the equations that describe the tidal 
evolution of these systems. 



2.1 Tidal evolution 

We base ou r model of tidal evolution on that described 
in detail in iDobbs -Dixon. Lin fc Mardlingj l|2004l ). We re- 
produce the relevant parts here. If the spin axes of a star 
and planet are aligned, tidal energy will be dissipated at 
rates defined by the tidal quality factors and Q'p. Con- 
sider a planet of mass Mp and radius Rp orbiting a star 



of mass M* and radius _R«. If the system has a semima- 
jor axis a and angular frequency n = (G{M, + Mp)/a^)^^^ , 
the rate of change of the eccentricty of the orbit is given 
bv (lEggleton. K iseleva fc Hut' 'l99d; iMardling fc LinI |2002| : 
iDobbs-Dixon. Lin fc Mardling.20oi ) 
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The stellar and planetary spin can also play an impor- 
tant role in the tidal evolution of close in planets. If the 
rotation axes are aligned and if M, >> Mp, the rate of 
change of the stellar spin frequency is (|Mardling fc Linll2002l : 
iDobbs-Dixon. Lin fc Mard ling 20di) 
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and cJ, is the change of stellar rotation due to angular mo- 
mentum loss through a stellar wind, which we will discuss 
in more detail in a later section. The quantity e, is the frac- 
tion of the mass of the star participating in tidal exchange, 
and Q, is the star's moment of inertia coefficient. We use 
e, = 0.1, appropriate for G stars which have shallow con- 
vection zones, and a* = 0.1. 

The rate of change of the planet's spin is 

2 \/,r\/r,\3 



Rr. 



/3(e) -/4(e) 



The quantities ep and Op are the fraction of the planet's 
mass involved in tidal exchange, and the planet's moment 
of inertia coefficient. Since gas giant planets can be fully 
convective, we use = 1 and use ap — 0.2. 

The total angular momentum perpendicular to the orbit 

is 

Jo{a,e) + Jpi^p) + J,{n,). (9) 

The orbital angular momentum is Jo{a,e) = MpMta?n{l — 
e^Y^"^ /{Mp + M,), while the angular momentum of the star 
and planet are 



(8) 



J, = 



ase^M^Rtn* 

2f 



Jfp — c^p J\^p Rp , 



(10) 

where it is assumed that only a fraction e, of the star, but all 
of the planet, is involed in angular momentum exchange. The 
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rate of change of the total angular momentum is then equal 
to the rate at which the system loses angular momentum 
through the stellar wind, Jui, ■ Differentiating equation 
and using J^^, — Q:*e,M,i?^iJ, gives 



Q,e.A/*i?*aJ. = Jo ( — 



2a 1 - 



+ 



(11) 



The rate of change of the semi-major axis of the orbit can 
then be determined by simply rearranging equation 



2.2 Stellar wind model 

We take the ste llar wind model from 

ICoUier Cameron fc JianQ l|l994 ). In the unsaturated 
dynamo regime with purely thermal driving, the rate of 
change of a star's angular frequency is 



where 
2 
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For the solar values we assume a magnetic flux density of 
Bo,0 = 3 G, a convective turnover time of Tc,q = 8.9 x 10^ 
s, and an angular velocity of ilo,0 = 4.0 x 10® rad s^^. To 
dete rmine Tc we assume a linear fit t o the values in Table 
1 of ICollier Cameron fc Jiankd l|l994h . The quantities rrip 
and fee are the proton mass and Boltzmann's constant and 
from lMestel fc SpruitI (|l987l ') we assume /? = 0.16. The final 
quantity k is the effective radius of gyration of the radiative 
interior and co nvective envelope combined, wh ich we take 
to be k^ = 0.1 (|Collier Cameron fc Jia"nk3ll994h . 
In the case of stars with saturated dynamos, 



(14) 



where Q. is the saturation limit w hich we get by assuming lin- 
ear fit to the values in Table 5 of lCollier Cameron fc Jiankd 
l|l994h . 



2.3 Basic setup 

We carry out Monte Carlo simulations of a large sample 
of close-in planets (a < 0.3 au) around solar-like stars. We 
consider a range of different initial radial distributions and 
eccentricity distributions that we will discuss in more de- 
tail in a later section. There are, however, some common 
elements to all of the simulations that we discuss here. 

Based on the properties of the known e xoplanets we 
distribute the planet masses as Mp oc M~^'^^ (|Marcv et al.l 
l2008 h with a minimum planet mass of Mp = 0.1 Mjup and 
a maximum planet mass of Mp = 15 Mjup. The minimum 
and maximum masses are chosen such as to represent the 
rang e of planet masses that cou ld undergo Type II migra- 
tion ([D'Angelo. Kiev fc Hennine 2003). The planet radii are 
computed using R oc M^-™ where m is varied from m = 
for low-mass planets to m = 1.5 for the most massive gas 
giants, to give a mas s-radius relation that matches that ex- 
pected theoretically ijChabrier et al.|[2009l l and that is nor- 
malised to give Rp = l Rjupiter when Mp = 1 Mjupitcr- The 
stars are assumed to have masses evenly distributed between 



0.75 M© and 1.25 M0 and have rotation periods distributed 
as a Gaussian with a peak at 8 days and with a half-width 
of 1.5 days. The planets are all assumed to start with rota- 
tion rates 10 times greater than their orbital periods (i.e., 
VLp = lOn). 

In each simulation we consider 1800 planets. The initial 
semi-major axis, eccentricity, planet mass, stellar mass, and 
stellar rotation period are all chosen randomly, but in such a 
way that the final distributions match those described above 
(we will discuss the chosen semi-major axis and eccentric- 
ity distributions below). The star and planet radii are then 
determined as above and the star mass and rotation peri- 
ods determine - using equations (|12p and (|13|l - the rate 
at which the star loses angular momentum through a stel- 
lar wind. For the star and planet, we use the tidal quality 
factors of Q', = l(f and = 5 x lO''. These are similar 
to those used by Jackson. Barnes fc Greenbergj l|2009i l and 
iDobbs-Dixon. Lin fc Mardlind (|2004l) a nd consistent with 
those determined by Brown et al.l I 201l|). We d id try using 
the tidal quality factors estimated by HansenI (2010), but 
these resulted in circularisation timescale for planets with 4 
and 5 day periods that were much longer than expected. 

Equations ([T]), ((Sj, ((8)1, and (pTjl - with equation 
pi|l rewritten in terms of d - are then integrated us- 
ing a fourth-order Runge Kutta integrator for a time 
chosen randomly between 2 x 10^ years and 6 x 10^ 
years. Any planet that reaches its Roche limit, given 
by a R = \Rp/0AQ2){M^/Mpf'^ (|Faber. Rasio fc Willemd 
l2005h. is assumed to be tid ally disrupted and destroyed 
I Jackson. Barnes fc Greenbera.200a '). 



3 RESULTS 

We assume here that close-in giant planets are unable 
to form in situ, but rather form beyond the snowline 
(a ^ 2.7 au for a Solar- mass star) and migrate inwards 
to th e ir current locations (|Lin. Bodenheimer fc RichardsonI 
1 19961 : iTrilling et all 1 19981 ). Although there are a number 
of mechanism that could result in planet migration, in 
the case of giant planets the most likely mechanism is 
Type II migration in which the p lanet migrates within 
a ga p in the proplanetary disc (|Goldreich fc Tremaind 



Il980l : iLin fc Papaloizoul 1 19861). Using theoretical models 



of Type II mig ration lArmitagd (|2007l ) (following on from 
lArmitage et~an (|2002l )) has shown that the resulting radial 
distribution of extrasolar planets is consistent with that 
observed, the best fit occuring if migration is assumed to 
be slightly suppressed at small radii. Although this ignores 
some effects, in particular planet-planet scattering which 
can also change the orbit al el ements of surviving planets 
(|Ford. Havlickova fc Rasiolf200ll ). it is at least a reasonable 

starting point. 

Inside 1 au, the best-fit model of lArmitagd (|2007l ) has 
a radial distrbution of dN/dloga oc a" *. We therefore use 
this, shown in Figure [T] as the base distribution for all of 
our simulations. We consider a few different initial eccen- 
tricity distributions, but since the final radial distributions 
don't depend strongly on the initial e ccentricity distribution 
jjackson. Barnes fc Greenberd 12009) , most of our simula- 
tions use the intial eccentricity distribution shown in Figure 

H 
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Figure 1. Initial radial distribution determined by assuming 
Type II in ward migration with migration slightly suppressed at 
small radii l lArmitagelbOOTl'l . Also shown is the curve dN/d log a oc 
a'' '* for comparison. 




e 



Figure 2. Initial eccentricity distribution used in the simulations 
presented. The final radial distributions do not, however, depend 
particularly strongly on the initial eccentricity distribution. 

3.1 Type II migration only 

Our first set of simulations consider the case o f Type II 
migra tion with slight suppression at small radii (|Armitagel 
We therefore use the initial distribution shown in 
Figure [T] and integrate each system for a randomly chosen 
time of between 2 and 6 Gyr. The final radial distribution 
is shown in Figure [S] This illustrates that the tidal inter- 
action causes most of the planets inside 0.04 au to migrate 
inwards ultim ately overfl owing their Roche lobes and be- 
ing destroyed (jjackson. Barnes fc Greenberg|l2009l ') . Outside 
0.04 au, however, the distribution is not significantly differ- 
ent to the initial distribution. 

The dashed-line in Figure [3] shows the radial distribu- 
tion of known close-in exoplanets. Since we are only inter- 
ested in those that are likely to undergo Type II migra- 
tion, we only include exoplanets for which M ^ 0.1 Jupiter 
masses. We have included only those plan ets discovered us- 
ing the doppler radial velocity technique (*Mavor fc Quelo3 
ll995l : IUdrv. Fischer i^Q ucloz 20o3). We considered also in- 
cluding planets detected via transits, but were concerned 
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Figure 3. The final semi-major axis distribution for the simu- 
lated planets (solid line) compared to the semi-major axis dis- 
tributed for observed exoplanets (dashed line) in systems with 
properties comparable to those considered in the simulations. 
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Figure 4. Cumulative semi-major axis distribution for the sim- 
ulated exoplanets (solid line) and for the observed exoplanets 
(dashed line), showing that the radial distribution of the observed 
exoplanets does not appear to be consistent with the simulated 
systems that are taken to have an initial radial distribution re- 
sulting from Type II migration alone. 



that selection effects may somewhat enhance any pile-up at 
small radii. So that our simulated sample and the observed 
sample can be compared, we have also excluded any planets 
in our simulated sample that would induce a stellar radial ve- 
locity of less than 2 m s~^, and only included those observed 
exoplanets whose host star masses are between 0.75 and 1.25 
Mq. Figure [3] suggests that these two distributions are not 
the same, with some evidence for a pile-up of observed exo- 
planets at around 0.05 au, that is not seen in the simulated 
distribution. Figure |4] shows the cumulative distribution for 
the two samples and again illustrates that a larger fraction of 
the observed exoplanets are located inside 0.1 au, compared 
to the simulated systems. Formally a Kolmogorov-Smirnoff 
test shows that the probability that the two distributions 
are drawn from the same parent distribution is Pks ~ 0.36. 
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Figure 5. The initial radial distribution of planets assuming 
Type II inward migration followed by a pile-up of planets inside 
0.1 au with the peak of the pile-up occuring where the orbital 
period is half the stellar spin period (a ~ 0.05 au), and a width 
corresponding to a variation in the orbital period of ±2 days. The 
three cases shown are one in which the number of planets inside 
0.1 is the same as in Figure[T] (solid line), an enhancement of 20 % 
over that in Figure [T] (dash-dot line), and a further enhancement 
of 30 % (dotted hue). 



3.2 Type II migration with stopping mechanism 

It hcis been suggested that rather than migrating all the 
way in to the central star, planets may be prevented from 
doing so by some kind of stopping mechanism that halts 
planets with periods of 4 days, corresponding to a semi- 
major axis of 0.05 au. One possibility is that the mag- 
netic fields of TTauri stars may be strong ei rough to dis- 
rupt the inner edg e of protoplanetary discs (|Konig]| [l99ll : 
iBouvier et al.ll2007 ). and will truncate the disc close to the 
corotation radius. Planets will continue to migrate inside 
this cavity, but this is expected to slow and become very 
inefficient once planets pa ss inside the 2:1 resonance with 
the inner edge of the disc (iLin. Bodenh eimer fc RichardsonI 



Figure 6. The final radial distribution of the simulated systems 
with initial radial distributions as shown in Figure[5] (solid, dash- 
dot and dotted lines ) compared to the observed radial distribu- 
tion of exoplanets (thick dashed line). 



20081) ■ The 



1 19961 : iKuchner fc Lecaij|2002l : iRice fc Armitag ^ 
distribution of stellar rotation periods (Hcrbst 'et al.ll2007h . 
with a peak at 8 days, would therefore produce a pile- up 
of planets with orbital periods ~ 4 days, corresponding to 
an orbital radius of ~ 0.05 au around a sola r-like star. 

As discussed above, lArmitaed (|2007l ) suggests that 
Type II migration would lead to a radial distribution, within 
1 au, approximated by dN / dloga oc aP'"^. This distribution 
would suggest that ~ 55 % of the planets between r — 0.01 
and r = 0.3 au would be located insided 0.1 au. We therefore 
consider 3 different pile- up scenarios, illustrated in Figure [5] 
The solid fine is a case in which the number of planets be- 
tween r = 0.01 and r = 0.1 au is the same as would be 
expec ted from the radial distribution derived by lArmitag^ 
l|2007h . The dash-dot line increases the fraction to 65 %, and 
the dash-dot-dot line is a case in which the fraction inside 
0.1 au is enhanced to ~ 88 %. 

As before, we evolve each system for a randomly cho- 
sen time of between 2 and 6 Gyr. The final radial distribu- 
tions are shown in Figure [S] The thick dashed- line is again 
the radial distribution of the known exoplanets and again 
we have only included those discovered using the doppler 



radial velocity technique and that orbit stars with masses 
between 0.75 and 1.25 Mq. We have also removed any simu- 
lated planet that would induce a radial velocity of less than 
2 m s"'^. The line styles in Figure [§] correspond to those 
shown in Figure [5] In all three case, the pile up remains 
but the best fit with the observed distribution is somewhere 
between the solid and dash-dot line. The cumulative distri- 
butions are shown in Figure [7] The line styles again corre- 
spond to those in Figure [5] and Figure |6] Again this illus- 
trates that the best fit to the observed distribution would 
be somewhere between the solid line and the dash-dot line. 
The Kolmogorov-Smirnoff test shows that both the solid and 
dash-dot lines, when compared with the dashed line, have 
a probability of Pks = 0.94 of being drawn from the same 
parent distribution. 

The above result suggests that, to explain the current 
distribution of close-in giant planets, there must have been 
a primordial pile-up of giant planets with the peak of the 
pile-up occuring at ~ 0.05 au. The best fit also occurs if 
the number of planets between 0.01 and 0.1 is similar to, 
but slightly higher than, that expected based on the ra- 
dial distribution of planets resulting from Type II migration 
(| Armitagel l2007h . One has to be slightly careful with this 
comparison, how ever, as the radial profile determined by 
lArmitaed (|2007| ) is a steady-state profile and does not indi- 
cate what fraction of planets will have migrated in side 0.1 au 
and then been lost. Population synthesis models (|lda fc LinI 
suggest that a large fraction of giant planets should 
migrate close to their host stars. If the stopping mechanism 
were efficient, we would then expect a significant primordial 
pile-up at small radii. Figures [6] and [7] are consistent with a 
pile-up but are not consistent with a particularly significant 
pile-up. This suggests either that the stopping mechanism 
is not efficient and that planets are lost before tidal evolu- 
tion becomes important, or that a smaller fraction of giant 
planets migrate close to their parent star than predicited by 
population synthesis models. 
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Figure 7. Cumulative semi-major axis distribution for the simu- 
lated systems (solid, dash-dot and dotted lines) with with initial 
radial distributions shown in Figure [S] compared to the cumula- 
tive radial distribution of the observed exoplanets (thick dashed 
line). This appears to illustrate that the observed distribution is 
consistent with a primordial pile-up at r ~ 0.05 au but not if there 
is a significant enhancement in the fraction of planets inside 0.1 
au. 
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Figure 9. The final radial distribution of the simulated systems 
with an initial radial distributions as shown in Figure \5\ (solid, 
dash-dot and dotted lines ) and with masses Mp > 2 Mjyp, com- 
pared to the observed radial distribution of exoplanets with the 
same range of masses (thick dashed line). The simulated systems 
still show a peak at ~ 0.05 au that isn't present in the observed 
systems. 
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Figure 8. The final radial distribution of the simulated systems 
with initial radial distributions as shown in Figure[5] (solid, dash- 
dot and dotted lines ) and with masses Mp < 2 Mjup, compared 
to the observed radial distribution of exoplanets with the same 
range of masses (thick dashed line). 



3.3 Mass and eccentricity dependence 

In Figure[S]we consider the entire mass range from Mp = 0.1 
Mjup to a maximum of Mp = 15 Mjup. If, however, we divide 
this into low-mass regime {Mp < 2 Mjup) and a high- mass 
regime (Mp > 2 Mjup) we find (Figure [S]) that for the low- 
mass planets the peak at a ~ 0.05 au remains and is still well 
matched for our simulated populations that have modest 
pile-ups at this radius. For the high-mass planets (Figure 
this is no longer the case. In the observed population there 

are no planets with a ~ 0. 05 au. 

It has been suggested l|Rice &: Armitag3l2008l ) that the 
interaction between a planet, located inside a magneto- 
spheric cavity, and the surrounding disc can lead to sub- 
stantial eccentricity growth for planets with AI > 1 Mjup 



and that the growth rate is greate st for the highest mass 
planets. It was therefore suggested (|Rice fc Armitag3l2008l ) 
that this would result in higher-mass, close-in planets being 
preferentially destroyed when compared with lower-mass, 
close-in planets. We therefore consider a situation in which 
the initial radial distribution is as shown in Figure [S] but 
the initial eccentricities increase with increasing mass. We 
assume, somewhat arbitrarily, the planets with masses be- 
low 1 Afjup have an eccentricity peak close to e = 0, but 
have a tail that can extend to e = 0.6, while planets with 
masses above lOMjup have eccentricities that lie preferen- 
tiall between e = 0.4 and e = 0.8. We assume a maximum 
eccentricity of e = 0.8 so that we don't simply lose a large 
fraction of the massive exoplanets immediately. We should 
note t hat simulations bv lBenftez-Llambv. Masset fc Beaugel 
l|201lf ) did not find eccentricity growth for massive planets 
in inner disc gaps. Their simulations, however, did not have 
co mpletely evacuated inne r cavities and so differ from those 
of iRice fc Armitagd l|2008l ). 

Figure \TU\ shows the resulting radial distribution for 
the lower-mass planets {Mp < 2 Mjup - solid line) which 
still agrees well with the observed radial distribution (thick 
dashed line). The radial distribution for the higher-mass 
planets {Mp > 2 Mjup - solid line) is shown in Figure [TT] and 
although not a particularly good fit to the observed distribu- 
tion (thick dashed line) does at least indicate that the pile- 
up is significantly depleted. Furthermore, if one considers 
the final eccentricity distribution, shown in Figure [131 again 
divided into planets with masses above 2 Mjupitor (solid line) 
and planets with masses below 2 Mjupitcr (dashed line), a 
group of high-mass exoplanets with large eccentricities re- 
mains and is at least qualitatively similar to that observed 
(Figure [12)) . This suggests that if higher-mass exoplanets do 
indeed preferentially have higher eccentricities than lower- 
mass exoplanets, this could explain why the pile-up at 0.05 
au is only observed for the lower-mass exoplanets and sug- 
gests that the mass dependent eccentricity distribution must 
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Figure 10. The final radial distribution of the simulated plan- 
ets with masses Mp < 2 Mjup and with an initial peaked radial 
distribution illustrated by the solid line in Figure [5] The initial 
eccentricity distribution is also assumed to depend on planet mass 
with the higher-mass planets preferentially having higher eccen- 
tricities than the lower-mass planets. For these lower- mass plan- 
ets the fit to the observed distribution (thick dashed line) is still 
good. 



1.2 : 

1.0 ; 
0.8 ; 

0.6 ; 

0.4 
0.2 



0.0 1 
0.01 



0.10 



Figure 11. The final radial distribution of the simulated plan- 
ets with mass Mp > 2 Mjup and with an initial peaked radial 
distributions illustrated by the solid line in Figure [5] The initial 
eccentricity distribution is also assumed to depend on planet mass 
with the higher-mass planets preferentially having higher eccen- 
tricities than the lower-mass planets. Although, for these higher- 
mass planets, the fit to the observed distribution (thick dashed 
line) is not particularly good, it is clear that the higher initial 
eccentricities has largely removed the initial peak at a ~ 0.05 au. 



be primordial rather t han be ing due to a difference in their 
subsequent evolution l|Rice fc Armitagai2008. V 



3.4 Planet-planet scattering 

The main goal of this work has been to show that an ini- 
tial distrib ution consistent with that expected from Type II 
migration l|Armitagell2007|) followed by a pile-up d ue to trun- 
cation of the inner disc jllice fc ArmitaedBoOSh is at least 
consistent - after tidal evolution - with the observed distri- 
bution of close-in exoplanets. We cannot, however, preclude 



Figure 12. Eccentricity distribution for high-mass (M > 2Mjup, 
solid line) and low-mass (M < 2Mj^p, dashed line) exoplanets, 
detected via radial velocity, with semimajor axes less than 0.3 au. 
There is a suggestion of a group of higher-mass exoplanets that 
retain relatively large eccentricities. 
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Figure 13. Simulated final eccentricity distribution for low-mass 
(M < 2Mjup, dashed line) and high- mass {M > 2Mjup, solid 
line) planets, with an initial mass dependent eccentricity distri- 
bution. In this case, a group of high-mass planets with large ec- 
centricities remains and is at least qualitatively consistent with 
that observed. 



the possibility that migration occurs through other mecha- 

nisrnS; 

IWright et a.l (|2009l ) have shown that there is a signifi- 
cant difference in the semi-major axis distribution of single 
and multiple-planet systems. The single-planet systems have 
a pile-up at ~ 3 days and a radial distribution that increases 
with increasing radius. The multiple-planet systems, on the 
other hand, have a semi-major axis distribution that is only 
weakly dependen t on r adius. It has been suggested (e.g., 
iMatsumura et al] l|2010l )) that this difference could be due 
to single-planet systems being dominated by planet-planet 
scattering (with additional planets being ejected from the 
system) while multiple-planet systems are dominated by mi- 
gration in a gas disc. Type II migration would, according to 
lArmitagg (200'f), result in a semi-major axis distribution 
that increases significantly with increasing radius. Qualita- 
tively, this is more consistent with the radial distribution of 
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the single-planet systems than it is with the multiple-planet 
systems, which have a semi-major axis distribution that is 
only weakly dependent on radius. 

It is hard to predict what kind of semi-major axis and 
eccentricity distribution planet-planet scattering would pro- 
duce, but we assume that if the distribution of single-planet 
systems is dominated by planet-planet scattering, and if the 
additional planets are ejected or are now beyond ~ 5 au, the 
eccentricity of the close-in planets must be large. We con- 
sider an initial eccentricity distribution that peaks at e 0.6 
and that has a half-width of 0.1. We consider 3 different ini- 
tial radial distributions, that shown in Figure [T] and 2 from 
Figure [5] (solid line and dotted line). We integrate these sys- 
tems in the same way as before. The results are shown in 
Figure 1141 The initially unpeaked distribution is shown as 
the solid-line. The initially weakly peaked distribution is the 
dotted line, while the initially strongly peaked distribution 
is the dash-dot line. A larger fraction of the close-in plan- 
ets are destroyed when compared with the simulations in 
which the initial eccentricity peaked at e = and so the 
initially weakly peaked distribution doesn't show a peak af- 
ter tidal evolution. That the most peaked initial distribution 
(dash-dot line) has a pile-up that is similar to that observed 
at least suggests that there could be an initial distribution 
that could result in a match with the observed distribution. 
However, tidal evolution also broadens the peak which sug- 
gests that the initial peak would need to be narrower than 
we assume. 

Although the above results do not preclude the 
possibility that the radial distribution of the observed 
single-planets systems which includes a peak at ~ 3 
days is due primarily to planet-planet scattering, it 
does suggest that this would require that planet-planet 
scattering is very effective at scattering planets into 
orbits with a < 0.1 au which is not entirely consistent 
with th eoretical expectations (iFord. Havlickova fc Rasid 



200ll: iMarzari F. fc Weidenschilhng 
Nagasawa, Ida fc Besshol i 20081 '). who 



20021 ). Even 
show that the 
Kozai mechanism can significantly enhance the number of 
close-in planets, conclude that main channel for forming 
close-in giant planets is probably still through Type II 
migration. Furthermore, in these simulations, the planets 
with a > 0.1 au undergo very little eccentricity evolution. 
That the observed sample of exopl anets has an eccent ricity 
distribution that peaks at e = (|Wright et a.l [20091 ) sug- 
gests that the initial eccentricity distribution, at least for 
those with a > 0.1, must have also peaked at e = which 
may suggest that the primary mechanism for getting these 
planets close to their parents stars cannot be planet-planet 
scattering unless something else (other than tidal evolution) 
then reduces their initial eccentricities. 



3.5 Mass-Period relation 

It has been suggested that the observed mass-period relation 
of close-in planets may g i ve us some information about their 
evolution. iFord fc Rasid (|2006l ) suggest that if planets evolve 
from very high initial eccentricities, their final orbital dis- 
tance will be abo ut twice the Roche distance (i.e., a ~ 2aij). 
iPont et al.l (|201ll ) suggest - by studying the mass-period re- 
lation of transiting planets - that tidal circularization and 
the stopping mechanism for close-in planets must be closely 
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Figure 14. The final radial distribution of the simulated sys- 
tems with initial eccentricity distributions that peak at e ~ 0.6 
and with initial radial distributions that are unpeaked (solid line) , 
weakly peaked (dotted line) and very strongly peaked (dash-dot 
line). Also shown is the observed radial distribution of exoplanets 
(thick dashed line). A peak only remains for the simulated sys- 
tem which was initially strongly peaked (dash-dot line) suggesting 
that, if the pile-up at ^ 3 days is primarily due to planet-planet 
scattering, it must be very effective at scattering planets into or- 
bits with a < 0.1 au. 



related. Our results here suggest that close-in planets are not 
typically circularised from highly eccentric orbits and also 
that any stopping mechanism occurs prior to tidal evolution 
playing a significant role. 

The initial semi-major axis distribution that, after tidal 
circularisation, produced the best fit when compared with 
the observed distribution of close-in planets was one in which 
there was a modest initial pile-up at a ~ 0.05 au (solid line 
in Figure[5|. In Figure[l5]we plot the planet-to-star mass ra- 
tio against orbital period for a randomly selected sample of 
these planets after they have undergone tidal evolution (as- 
terisks) together with the observed close-in planets detected 
via radial velocity (triangles). For the modelled systems we 
have randomly selected the same number of planets as the 
number of observed close-in planets. The two diagonal lines 
are the Roche limit and twice the Roche limit. Qualitatively, 
the mass-period relation of the modelled systems looks very 
similar to that of the observed systems. Even though our 
modelled sample was not tidally evolved from an initially 
highly eccentric orbit, the inner boundary is just beyond 
twice the Roche limit. Admittedly, the sample is quite small 
and we have not analysed this in extensive details. It does, 
however, appear as though evolution through Type II mi- 
gration (with a mechanism for stopping planets at a ~ 0.05 
au) followed by tidal evoluion can produce a mass-period 
relation that is consistent with that ob served. It has also 
been suggested l|Davis fc WheatlevllioogI ) that close-in plan- 
ets can also be lost through evaporation. Although evapora- 
tion may well be operating, the lack of planets close to their 
Roche limit is a natural consequence of inward migration 
followed by tidal evolution. 
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Figure 15. Scatter plot showing tlie mass-period relation for a 
randomly selected sample of the modelled systems (asterisks) to- 
gether with the mass-period relation for the observed systems 
(triangles). The model we have chosen is the one in which there 
was an initial modest pile-up at a ~ 0.05 au ( solid line in Figure 
\5^. For the observed systems, we have only considered those de- 
tected via the radial velocity technique. The two diagonal lines are 
the Roche limit and twice the Roche limit. The two populations 
appear, qualitatively, to be very similar. 



4 DISCUSSION AND CONCLUSIONS 

The primary goal of this paper was to determine if the cur- 
rent distribution of close-in, gas-giant {M > 0.1 Mjup) exo- 
planets around Solar-like (0.75 ^ M* 1.25 M©) stars 
is consistent with an initial distribution resulting primarily 
from gap-opening migration in a gas disc (Type II) (e.g., 
lArmitage etal] |2002| : lArmitagd |2007| ). The chosen initial 
distributions were evolved using tid al evolution equations 
(jPobbs-Dixon. Lin fc Mardlind 20041 ) for a randomly cho- 



sen time between 2 x 10^ and 6 x 10® years. 

The results suggests that an initial radial distribu- 
tion due to Type II migration alone does not match the 
currently observed radial distribution which shows a pile- 
up of planets at a ~ 0.05 au. If, however, the inner re- 
gions of the disc are truncated due to interactions with 
the star's magnetosphere, it has been suggested that plan- 
ets should pile- up when inside the 2:1 resonance with 
the inner disc edge (Lin, B odenheimer fc Richardsonlll996l : 
iRice fc Armit agc 2008) which would correspond to a semi- 
major axis of ~ 0.05 au for Solar-like stars with initial rota- 
tion periods of ~ 8 days. The present study shows that the 
mass-period-eccentricity relation is indeed consistent with 
expectations from disc migration with a magnetospheric gap 
and tides, at least for stars with masses in the range between 
0.75 Mq and 1.25 M© and as long as the primordial pile-up is 
slightly (but not significantly) enhanced compared to what 
would be expected based on a steady-state model of Type II 
migration. We also find a good, qualitative, agreement be- 
tween the mass-period relation of the modelled systems and 
that of the observed systems. 

This pile-up is, however, only evident for the lower- 
mass planets (Mp < 2 Mjup). Simulated systems with a 
mass-independent initial eccentricity distribution retain an 
initial pile-up for all planet masses. The pile-up for higher- 
mass planets is, however, significantly reduced if the ini- 



tial eccentricity distribution is assumed to be mass de- 
pendent with the higher-mass planets preferentially hav- 
ing higher eccentricities . This is consistent with simulations 
l|Rice fc Armitagdliooi ) suggesting that higher-mass plan- 
ets will have enhanced eccentricity growth inside a magne- 
tospheric cavity. 

We should stress that this does not prove that the pri- 
mary mechanism for producing the observed close-in giant 
exoplanets is Type II migration but simply shows that - for 
reasonable assumptions about the initial distributions - it 
is consistent with this being the case. We also consider one 
set of simulations with an initial eccentricity distribution 
that peaks at a large eccentricity, aimed to mimic a possible 
planet-planet scattering scenario. Although it can result in 
a final distribution that matches that observed, it requires a 
more significant initial pile-up when compared to an initial 
eccentricity distribution that peaks at e = 0. Furthermore, 
the planets beyond a ~ 0.1 au retain their high eccentricities 
which is not consistent with observations. 

The observation of misaligned , and in some 
cases retrograde, clo s e-in p lanets (IWiiin et al.l l2009l : 
ICoUier Cameron et al.l I2OI0I : iTriaud et al.l I2OI0I ) does, 
however, strongly suggest that planet-planet scatter- 
ing (which may undergo "Kozai" cycles if a binary 
or distant planetary companion is present) must play 
a role in the formation of some close-in exoplanets. 
It is intriguing that misaligned systems occur at all 
masses, but dom i nate for stars with Tbs > 625 K 
dWinn et al.l I2OIOI : iBarnes. Linscott fc Shporeij I2OIII ). In 
their po pulation-synthesis studies t hat in clude disc mi- 
grations lAlibert. Mordasini fc Bend (|201ll ) find that the 
short disc lifetimes of high-mass stars doesn't leave enough 
time for Jupiter-mass planets to migrate into close orbits. 
Our results, therefore, appear to be consistent with this. 
For stars with masses between 0.75 and 1.25 M0, the 
distribution of close-in giant planets can be modelled by 
assuming Type II migration followed by tidal evolution. 
There will be some systems that have been influenced 
by planet-planet scattering, which may have undergone 
"Kozai" cycles, but these do not dominate. For the more 
massive stars (A/* > 1.25 M0), where the disc lifetime is 
too short for Type II migration to be effective, close-in 
giant planets are preferentially formed through dynamical 
interactions in which "Kozai" cycles may also operate. 



ACKNOWLEDGEMENTS 

The authors would like to thank the anonymous ref- 
eree for helpful comments that have improved this paper. 
W.K.M.R. acknowledges support from the Scottish Univer- 
sities Physics Alliance (SUPA) and for support from the 
Science and Technology Facilities Council (STFC) through 
grant ST/H002380/1. The authors would also hke to thank 
the Isaac Newton Institute for Mathematical Sciences for 
their hospitality during the Dynamics of Discs and Planets 
Programme and would like to thank Phil Armitage for many 
useful discussions. 



REFERENCES 

Alibert Y., Mordasini C, Benz W., 2011, A&A, 526, A63 



© 0000 RAS, MNRAS 000, 000-000 



10 W.K.M. Rice, J. Veljanoski, & A. Collier Cameron 



Armitage P.J., Livio M., Lubow S.H., Pringle J.E., 2002, 

MNRAS, 334, 248 
Armitage P.J., 2007, ApJ, 665, 1381 

Barnes J.W., Linscott, E., Shporer A., 2011, ApJS, in press. 
Bell K.R.. Cassen P.M., Klahr H.H., Henning Th., 1997, ApJ, 

486, 372 

Bcnftcz-LLamby P., Masset F., Beauge C, 2011, A&A, 528, A2 
Bouvier J., Alencar S.H.P., Harries T.J., Johns-KruII CM., 

Romanova M.M., 2007, in Reipurth B., Jewitt D., Keil K., 

eds, Protostars and Planets V. Univ. Arizona Press, Tucson, 

p. 479 

Brown D.J. A., Collier Cameron A., Hall C, Hebb L., Smalley 

B., 2011, MNRAS, 415, 605 
Chabricr C, BarafEc I., Lccontc J., Gallardo J., Barman T., in 

Stempels E., ed., AIP Conf. Proc. Vol. 1094, Cool Stars, 

Stellar Systems and the Sun: Proe. 15th Cambridge 

Workshop on Cool Stars, Stellar Systems and the Sun. Am. 

Inst. Phys., New York, p. 102 
CoUier Cameron A., Jianke L., 1994, MNRAS, 269, 1099 
CoUier Cameron A., Bruce V.A., Miller G.R.M., Triaud 

A. H.M.J., Queloz D., 2010, MNRAS, 403, 151 
Davis T.A., Wheatley P.J., 2009, MNRAS, 396, 1012 
D'Angelo C, Klcy W., Henning T., 2003, ApJ, 586, 540 
Dobbs-Dixon I., Lin D.N.C., Mardling R.A., 2004, ApJ, 610, 464 
Eggleton P.P., Kiseleva L.G., Hut P., 1998, ApJ, 499, 853 
Faber J., Rasio F., Willems B., 2005, Icarus, 175, 248 
Fabrycky D., Tremaine S., 2007, ApJ, 669, 1298 

Ford E.B., Rasio F.A., 2006, ApJ, 638, L45 

Ford E.B., Havlickova M., Rasio F.A., 2001, Icarus, 150, 303 

Goldreich P., Tremaine S., 1980, ApJ, 241, 425 

Hansen B.M.S., 2010, 723, 285 

Herbst W., Eisloffel J., Mundt R., Scholz A., 2007, in Reipurth 

B. , Jewitt D., Keil K., eds, Protostars and Planets V. Univ. 
Arizona Press, Tucson, p. 297 

Ida S., Lin D.N.C., 2004, ApJ, 616, 567 
Jackson B., Barnes R., Greenberg R., 2009, ApJ, 698, 1357 
Konigl A., 1991, ApJ, 370, L39 
Kozai Y., 1962, AJ, 67, 591 
Kuchner M.J., Lecar M., 2002, ApJ, 574, L87 
Lin D.N.C., Papaloizou J., 1986, ApJ, 309, 846 
Lin D.N.C., Bodenheimer P., Richardson D.C., 1996, Nature, 
380, 606 

Marcy, G.W., et al., 2008, Phys. Scr., 130, 014001 
Mardling R.A., Lin D.N.C., 2002, ApJ, 573, 829 
Marzari F., Weidenschilling S.J., 2002, Icarus, 156, 570 
Matsumura S., Thommes E.W., Chatterjee S., Rasio F.A., 2010, 

ApJ, 714, 194 
Mayor M., Queloz D., 1995, Nature, 378, 355 
Mestel L., Spruit H.C., 1987, MNRAS, 226, 57 
Morton T.D., Johnson J.A., 2011, ApJ, 729, 138 
Nagasawa M., Ida S., Bessho T., 2008, ApJ, 678, 498 
Queloz D., et al., 2010, A&A, 517, LI 

Pont F., Husnoo N., Mazeh T., Fabrycky D., 2011, 414, 1278 
Rasio F., Tout C.A., Lubow S.H., Livio M., 1996, ApJ, 470, 1187 
Rice W.K.M., Armitage P.J., 2008, MNRAS, 384, 1242 
Triaud A.H.M.J., 2010, A&A, 524, A25 

Trilling D.E., Benz W., GuiUot T., Lunine J., Hubbard W.B., 

Burrows A., 1998, ApJ, 500, 428 
Udry S., Fischer D., Queloz D., 2007, in Reipurth B., Jewitt D., 

Keil K., eds, Protostars and Planets V. Univ. Arizona Press, 

Tucson, p. 685 
Ward W.R., 1997, Icarus, 126, 261 

Winn J.N., Fabrycky D., Albrecht S., Johnson J.A., 2010, ApJ, 

718, L145 
Winn J.N., et al., 2009, ApJ, 703, L99 

Wright J.T., Upadhyay S., Marcy G.W., Fischer D.A., Ford 

E.B., Asher Johnson J., 2009. ApJ, 693, 1084 
Wu Y., Murray N., 2003, ApJ, 589, 605 

© 0000 RAS, MNRAS 000, 000-000 



